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We have studied the analytical dynamics of Bell nonlocality as measured by CHSH inequality and 
entanglement as measured by concurrence for two noisy qubits that have dipole-dipole interaction. 
The nonlocal entanglement created by the dipole-dipole interaction is found to be protected from 
sudden death for certain initial states. 



I. INTRODUCTION 

Entanglement is a kind of quantum correlation and refers to separability of the states. The entangled states 
are the main resource for quantum information and computation applications, such as quantum teleportation [T], 
superdense coding [2J, and quantum cryptography [3]. Entanglement and other types of quantum correlations are 
not easy to maintain when the system is in contact with an environment which leads to decoherence (Jj. Thus the 
study of entangled states under decoherence is one the most important aspects in theoretical as well as experimental 

qv areas [5J [5] . It was shown that the interaction between the single qubit and its environment leads to exponential decay 

of the qubit coherences while entanglement between two such qubits can cease to exist in a finite time, a phenomenon 
which was named "entanglement sudden death" (ESD) by Yu and Eberly [7]. Later Lopez et al., revealed that 
when the bipartite entanglement suddenly disappears, the entanglement of the corresponding reservoir suddenly and 
necessarily appears, which is a counter phenomenon to ESD and termed "entanglement sudden birth" (ESB) [5]. For 
bipartite systems, many theoretical efforts have been devoted to ESD [7J I9HT4] as well as ESB [8j [15j [16] and also 
ESD is verified experimentally [6]. However, they need further clarifications for a deeper understanding. 

On the other hand, it was demonstrated by Werner that there are some entangled mixed states whose correlations 
can be reproduced by a classical local model [T7] . The only way of identification of such states is to use violations of 
Bell-inequalities. The violations of Bell inequalities identify the genuine multipartite entanglement which is useful for 
quantum computations [TSJ [I^]- For a bipartite system, there are many Bell type inequalities, such as quadratic- and 

£S1 CHSH-Bell inequalities [201 HI] • The CHSH-Bcll inequality is a good indicator of nonlocal correlations and the relation 

^- with the entanglement has been already known 22 29 . Among them, Bellomo et al., found that Bell inequality might 
not be violated for a state with high entanglement for a two qubit system subject to amplitude damping [22]. Kofman 
et al., showed that the survival time for entanglement should be much longer than the Bell inequality violation under 
dephasing and energy relaxation type channels at an arbitrary temperature |25) . And recently, Li et al., studied the 
violations of Bell inequality under classical Ornstein-Uhlenbeck noise and showed that strong non-Markovian effect 

C ■ can protect the nonlocal entanglement identified by the violation of the CHSH inequality [24] . 

In this paper, we have analyzed the dynamics of Bell nonlocality as measured by CHSH inequality and entanglement 
as measured by concurrence for a system of two coupled qubits that interact via dipole-dipole interaction. The effect 
of environment on the qubit is modelled as a stochastic energy level with Ornstein-Uhlenbeck type correlation [TU] . 
The effects of non-Markovianity of the dynamics, purity of the initial states as well as the dipole-dipole interaction 
strength on Bell nonlocality and entanglement have been investigated for a system initially prepared in extended 
Werner-like states. 

The organization of this paper is as follows. In Sec. [Hj we introduce the model and its solution by solving the 
master equation for the two qubit reduced density matrix. In Sec. |IH[ we briefly discuss the Wootters concurrence 
as well as CHSH Bell inequality and the effect of purity, dipole-dipole interaction strength and the non-Markovianity 
on entanglement and Bell nonlocality is studied for extended Werner-like initial states. In Sec. |IV[ we conclude as a 
summary of important results. 
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II. THE MODEL AND ITS SOLUTION 

In the present paper, we consider two interacting qubits that are coupled to their independent environments which 
leads to fluctuating energy levels. The qubit-qubit interaction is assumed to be of the type Heisenberg XX model. This 
model can be thought of as the Kubo- Anderson model extended to two coupled qubits [30] ■ The typical Hamiltonian 
for this model can be given as [PJHUJOII] (we set h— 1): 

tt r + \ j(aA*b j_'A^.B\ , Raw ^a , M*) ~b /-in 

Utot{t) = J{a x a x +v y Vy)+ 2 cr z H ^ — cr z , (1) 

where J is the qubit-qubit interaction strength, a i ' {i = x, y, z) are the usual Pauli spin operators and Qa,b (t) are the 
independent frequency fluctuations of the qubits with mean value properties that obey non-Markovian approximation: 

(2) 

■=A,B, (3) 

where M{...} stands for the statistical mean over the noise f^i(t) and f2s(t). Here Ti(i — A, B) are the damping rates 
due to the coupling to the environments, 7 is the noise bandwidth which determines the environment's finite correlation 
time (t c = 7 _1 ) and a(t — s) is the reservoir correlation function. For simplicity, we will take the noise properties to 
be the same for A and B (e.g., Ta = ^b = r). And note that in the limit 7^00 (r c — >• 0), Ornstein-Uhlenbeck noise 
reduces to the well-known Markovian case [TT] : 

a(t-s)=TS(t-s). (4) 

For the total system described by the Hamiltonian ([!]), the stochastic Schrodinger equation is given by 

if t mt)) =H tot (t)Mt)), (5) 

with formal solution 

\y(t)) = u(t,n A ,n B )\y(o)), (6) 

where the stochastic propagator U(t, £Ia,^b) is given by 

U(t,n A ,n B )=e-^oH t o t (s)ds^ ^ 

The reduced density matrix for qubits A and B is then obtained from the statistical mean 

= M{|tt(t))<tt(t)|}. (8) 

With the help of the raising and lowering operators, &±' = [p A,B ± ia A ' B )/2, and the stochastic Schrodinger 
equation ([5| , the master equation for the reduced density matrix can be derived as [J2 l3"2Tf3"4"] : 

j t p = -i[H, p] - ^(2/5 - & A pa A - afpaf), (9) 

where G(t) = L a(t — s)ds = ^(1 — e -7 *) is a time-dependent coefficient which includes the memory information 
of the environmental noise and H = 2J{a A <r B _ + a A a B ) is the interaction Hamiltonian which represents the dipole- 
dipole interaction between the qubit A and B [35] . Note that the first term on the right-hand side of Eq. ffity leads to 
oscillatory dynamics while the second term causes decay. 

The differential equations governing the time evolution of the system in the standard basis {|1) = |11) , |2) = 



|10) , 1 3) = |01) , 1 4) = 1 00)} can be easily calculated: 



Pii 


= (* = 1,4), 


P22 


= 2iJ(p 23 - p 23 ), 


P33 


= -2iJ(p23 - P23), 


Pl2 


= 2iJp 13 - G(t)pi2, 


Pl3 


= 2iJpi2 - G(t)pi3, 


pu 


= -2G(t)p 14l 


f>23 


= 2i,J( P22 - p 33 ) - 2G(t)p 23 


P24 


= -2iJp 3i - G(t)p 2 i, 


P3i 


= -2iJp 2i - G(t)pu, 



(10) 

where the asterisk in the superscript of p* mn denotes the complex conjugation of p m n- After a simple calculation, 
the analytical solutions of the coupled first order ordinary differential equations in Eq. ( 10 1 can be found. They are 
reported infAl It should be noted that in the limit of J — > 0, the solutions (Al) have simple analytical forms which 
are analyzed in Refs. [5] and [24"] . 



III. ENTANGLEMENT AND CHSH-BELL INEQUALITY DYNAMICS 



For two-qubit systems, Wootters concurrence can be used as a measure of entanglement [35]. The concurrence 
function varies from C — for a separable state to C = 1 for a maximally entangled state. The concurrence function 
may be calculated from the density matrix p for qubits A and B as: 

(11) 



C(p) = max{0, y/Xi - \J \ 2 - \J\ 3 - \J A 4 }, 
where the quantities A^ are the eigenvalues in decreasing order of the matrix 

Ptrans = K°y ® Vy)p*(°y ® df). (12) 

In the following, we consider entanglement dynamics of the qubits whose density matrix has a common X-form [llj : 



pu pu 

p 22 p 23 

p 32 p 33 

p 4 i pu 



(13) 



Such X-states arise in a wide variety of physical situations and include pure Bell states |37| as well as the well-known 
Werner mixed states [TTj- Note that the time evolution of Eq. ( 10 1 preserves the form of X-state. Then one can easily 
show that the concurrence function for the X-state ( 13 ) is given by [9] 



C(p) = 2max{0, |p 2 3| - \J PuPu, \pu\ - ^/P22P33~}- 
For two-qubit case the CHSH-type Bell inequality can be written in the following form 

B(p) - \{O a O b ) - (OaO' b )\ + (0' A B ) + (0' A 0' B ) < 2, 



(14) 



(15) 



where (OaOb) = Tr(pOAOs) is the correlation function and Os = Os.cs, wherein Os = 

(sin 6>5 cos 05, sin 6*5 sin 0s, cos 6*5) is the unit vector, as — (erf, erf, erf) is the Pauli matrices vector and 0' s — 

Os{^'s7 ( f ) s)- Bcllomo et al, showed that the maximum of the Bell function B(p) for a X-structured density ma- 
trix (|13l) can be given as [HI [53] : 



B, ; 



,(p) = 2^pi(t) + Q2(t), 



where 



P(t) 

Q(t) 



Pu + Pu - P22 
2(| /0l4 | + |p 23 |)- 



P33, 



(16) 



(17) 



It should be noted that in the regions where B maK (p) > 2, the CHSH inequality is violated. It implies that the state 
p is genuinely bipartite Bell nonlocal, thus the correlations cannot be accessible by any classical local model. Also 
notice that the expression for -B max (/5) in Eq. (16) coincides with the one that would be obtained using the formal 
Horodecki criterion [55] , 

In the following, we consider two classes of initial states, called extended Werner-like (EWL) states in the form: 
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(18) 



+rry 



where r is the purity of the initial states which ranges from for maximally mixed states to 1 for pure states, 1^ is 
the 4x4 identity matrix, |$) = a 1 10) + f3 |01) , |\&) = a |11) + /3 |00) are the Bell-like pure states and the parameter 
a is sometimes called the degree of entanglement (here we set a as real number and f3 = y/1 — a 2 ) |22j . From 
Eq. (10) or (Al ), it can be easily noted that the initial states of Eq. (fl8|) belong to "X" states at any time t. So the 



time-dependent concurrence and the maximum of the Bell function can be obtained from the the Eqs. (14) and (16), 
respectively. 

Before starting our qualitative analysis, we want to emphasize that the time-dependent concurrence and the max- 



imum of the Bell function for /3*(0) have simple analytical forms. According to Eq. (Al), for this initial state the 
density matrix at any time can be found as: 
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(19) 



then the time-dependent concurrence and the maximum of Bell function for the p^{0) initial state can be obtained 



as: 



C*(p) = 2max{0,ra/3e- 



■2/(t) 



b: 



,05) = 2ry / l+4(a/3e- 2 /( t )) 2 , 



(20) 



where /(£) = | it + Me 7t — 1) 1 . For the initial state /3$(0), a similar expression to Eq. (20) can be written by using 

the solution in [A] But the elements are quite involved and we do not display them here for brevity; p<s>(t) for special 
values of a and J are displayed and discussed below. One should note that the dynamics of p*(0) initial state is 
undisturbed by the dipole-dipole interaction; which is expected because of the single excitatio n na ture of dipole-dipolc 



interaction and the form of |\1/) state which is |\I>) = a |11) + yl — a2 1 00) . Furthermore, Eq. (20) also corresponds to 
the time-dependent analytical expressions of C(p) and B max (p) for the initial states /5$(0) and p*(0) in the absence 
of dipole-dipole interaction (i.e., J = 0) |24|. The dynamics of entanglement and Bell nonlocality for the initial state, 
(5^(0), was extensively examined in Ref. [24], thus we will not cover this state in detail at the rest of this paper. 

The dynamics of concurrence and Bell nonlocality as a function of the dipole-dipole interaction strength and the 
dimensionless time are displayed in Fig. 1(a) and (b) for a = 1 and in Fig. 1(c) for a = 1/V4 for the initial states 
p$(0) with r = 0.95. All figures and quantities show a faster than exponentially decaying oscillatory behaviour with 
J-dependent oscillation frequency. Based on Eq. (20), for a = 1, it should be noted that in the absence of dipole-dipolc 
interaction (i.e., J = 0), the states p^{t) and ^(tjhave no entanglement and do not violate CHSH inequality at any 
time, while in the presence of dipole-dipole interaction (i.e., J 7^ 0), p$(i) possesses a high degree of entanglement and 
violates the CHSH-Bell inequality as can be seen from Fig. 1(a) and (b). On the other hand, for a = l/-\/4, the initial 
state (5$ (0) is entangled at t — and the entanglement as well as Bell nonlocality for J ^ live longer compared to 
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FIG. 1: Fig. 1(a) and 1(b) show C(p) (solid plot) and B mayc (p) - 2 (dashed plot) versus Ft for a = 1 and J = 0.4 (Fig. (a)) 
and J = 0.8 (Fig. (b)). Fig. 1(c) shows C(p) (thick plot) and B max (p) — 2 (thin plot) versus Tt for a = l/\/4 and J = (solid 
plot), J = 0.4 (dashed plot) and J = 0.8 (dotted plot). Here the figures are plotted for p$(0) with r — 0.95 and 7/T = 0.1. 

the J = case (see Fig. 1(c)). The common feature of all figures is that Bell nonlocality is found to die out before 
the entanglement does as also found by Kofman and Korotkov in the case of dephasing and dissipative dynamics of 
two qubits [25) . 




FIG. 2: C(p) (Fig. (a)) and B max (/5) - 2 (Fig. (b)) versus Ft and -y/T for p&(0) with J = 0.5, a = 1 and r = 0.95. 

Non-Markovianity dependence of the entanglement and Bell inequality violation for /5$(0) with a = 1 is displayed 
in Fig. 2(a) and 2(b), respectively. For the relatively high purity (r = 0.95) considered in this case, the effect of non- 
Markovianity is found to be a prolonging of the lifetime of both concurrence and the Bell inequality violation. This 
finding is in line with many similar studies on the effect of non-Markovianity on the dynamics of entanglement, since 
non-Markovianity is a measure of memory effects, its increase may lead to a longer lived or protected entanglement [51 
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FIG. 3: C(p) (Fig. (a) and (b)) and B max (p) - 2 (Fig. (c) and (d)) versus Ft and r for p$(0) with J = 0.5, a = 1. Here Fig. (a) 
and (c) correspond to Markovian dynamics with 7/T = 10 and Fig. (b) and (d) to non-Markovian dynamics with 7/T = 0.1. 



The effect of purity of the initial state on the dynamics of p$ (£) state for Markovian and non-Markovian environments 
is displayed in Figs. 3(a)-(d) for the concurrence and Bell inequality with a = 1. It should be noted that the dipole- 
dipole interaction creates a significantly wide purity range for non-zero entanglement as 0.4 < r < 1 for Markovian and 
0.35 < r < 1 for non-Markovian dynamics. On the other hand, the purity dependence of non-zero Bell nonlocality is 
narrower compared to the entanglement as 0.8 < r < 1 and 0.7 < r < 1 for Markovian and non-Markovian dynamics, 
respectively. The figures also show a pronounced difference between the dynamical behaviour of the both quantities 
for the pure and the mixed states. While both entanglement and Bell nonlocality decreases exponentially with time for 
pure states under Markovian and non-Markovian dynamics, there is a sudden death phenomenon for both quantities 
for the initially mixed states. 

For the special case of a = l/\/2, the initial states considered here become 



P*(0) = ^/ 4 + r|$)($|, 
M0) = ^h + r\*){*\, 



(21) 



where |$) = -4(|10) + |01)) and |#) 
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presented in g the time-dependent density V matrix for 0.(0) can be written as 



|11) + 1 00)) are the pure Bell states. With the help of the solutions 
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(22) 



then the concurrence and the maximum of the Bell function for the p$(t) of Eq. (22 1 can be expressed in a simple 



analytic form as 



B, 



C*(p) = max{0,re- 2/( * 
(p) = =2r\[\- 



4/(*) 



(23) 



For the special value of a = l/v2, the dynamics of the entanglement and the Bell nonlocality were investigated by 
Li and Liang who showed that C(p) and B max (p) show exactly same dynamics for p$(0) and p^ (0) initial states in 
the absence of dipole-dipole interaction [23] , thus we will not discuss it here. 




FIG. 4: C(p) (Fig. (a) and (b)) and B max (/3) - 2 (Fig. (c) and (d)) versus Vt and a for p # (0) or p*(0) with J = 0,7/T = 0.1 
and r = 0.95 (Fig. (a) and (c)) and r = 1 (Fig. (b) and (d)). 

To further elucidate the role of the so called degree of entanglement parameter a in the dynamics of C(p) and 
-Bmax(/5); w e consider the general initial states given by Eq. (18). It is easy to write down the time-dependent density 
matrix elements for the type of initial states of Eq. (18) by using the solution in [A] The dynamics of concurrence 
and Bell nonlocality as a function of the parameter a of the initial state are presented in Fig. 4 and 5. In Fig. 4, 
we display the behavior of C{p) and B mstx (p) — 2 for the non-interacting case for the initial state /5$(0) or p^ (0) , 
while Fig. 5 presents the dynamics for the J = 0.5 value of the dipole-dipole interaction strength for /5$(0). For the 
non- interacting case, the effect of purity and the a parameter on the dynamics of C{p) and -B m ax(/S) — 2 are displayed 
in Figs. 4(a) and (b) and Figs. 4(c) and (d), respectively. As can be seen from these figures, the entanglement and 
Bell nonlocality are maximum for a — \j\f2 which corresponds to initial Bell-state which offers the longest time of 
survival for the mixed state case (Fig. 4(a) and (c)) in the absence of dipole-dipole interaction. The volume of nonzero 
entanglement and Bell nonlocality decreases as r (purity) decreases as can be seen from a comparison of Figs. 4(a) 
and 4(b) and also 4(c) and 4(d). For r = 1, the entanglement and Bell nonlocality have non-zero values for all a values 
with exponential decay, except for a = or a = 1. One should also note that Figs. 4(a)-(d) are for non-Markovian 
dynamics, in the Markovian case Fig. 4(b) and (d) would be qualitatively similar, but the nonzero regions in Fig. 4(a) 
and (c) would shrink EH H3]- 

The effect of dipole-dipole interaction on the dynamics C{p) and B mlix (p) — 2 as function of a and Tt are displayed 
in Figs. 5(a)-(d) for p$(0) with J = 0.5. Comparing to the corresponding subfigures of Fig. 4 and 5, the most 
pronounced difference is in the dynamics of the initial state with a = I/a/2 which suffers the quickest death among 
all the other probable initial mixed states. Also a = and a = 1 states which are initially non-entangled are found 
to have entanglement for the longest time under dipole-dipole interaction. Also we should mention that for the pure 
initial states (r = 1), the dipole-dipole interaction leads to non-zero entanglement and violation of Bell inequality for 
a = and a = 1 for all times as can be seen from Figs. 5(b) and 5(d). 
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FIG. 5: C(p) (Fig. (a) and (b)) and B max (p) - 2 (Fig. (c) and (d)) versus Yt and q for p # (0) with J = 0.5, 7/r = 0.1 and 
r = 0.95 (Fig. (a) and (c)) and r = 1 (Fig. (b) and (d)). 



IV. CONCLUSION 



We have investigated the dynamics of entanglement as measured by concurrence and Bell nonlocality as measured by 
CHSH inequality for two noisy qubits which are connected to each other by dipole-dipole interaction. The considered 
noise have Ornstein-Uhlenbeck type correlation and the dynamics is investigated analytically for extended Werner- 
like initial states. We have found that the most important effects of dipole-dipole interaction to the dynamics of 
the considered quantities are as follows: The dipole-dipole interaction can create nonzero entanglement and lead to 
violation of Bell inequality for /5$(0) = -^p/ 4 + r |$) ($| initial state (where |$) = a |10)+\/l - a 2 |01) ) with a = or 
a = 1 which can be protected from sudden death for a pure initial state (r = 1). On the other hand, for mixed initial 
states (r < 1) the dynamics of p$(0) with a = l/v2 is found to have longest life time for nonzero entanglement and 
Bell nonlocality violation in the absence of dipole-dipole interaction. The considered quantities suffer the quickest 
death for a = l/v2 if the dipole-dipole interaction between qubits is considered. 



Appendix A 

The solutions of the differential equations (flOl) can be obtained as 



where 



pa = Pii(0) (i = 1,4), 

Y + Z 
P22 — — g — ' 

Y - Z 

P33 — 2 ' 

P12 = (pi 2 (0)cos(2Ji) + zp 13 (0)sin(2Ji))e- / W, 

Pis = (pi 3 (0)cos(2Ji) + zp 12 (0)sin(2Ji))e- / ^, 

Pl4 = Pu(0)e- 2/ W, 
A + B 

P23 - -^— , 

P24 = (p 24 (0)cos(2Jt)- ip 3i (0)sm(2 Jt))e~ m , 

P34 = ( / 9 3 4(0)cos(2Jt)-ip24(0)sin(2Jt))e- / W, (Al) 
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A = (p 2 3(0)+p; 3 (0))e- 2 ^, 
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V 7 7 
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V 7 7/ 
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V 7 7/ 
+ 8iJ(p 23 (0)-p 3 2(0))), 

i \" + _ / . e r\ ,„ , 2 2x „ / _ e r 
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V 7 7/ V 7 7 

x ((p 32 (0) - p 23 (0))(r + e) - 8zJ(p 22 (0) - ^(O)))), 

£> = (^J iFi(«+;i + -i — ^J( 7 r(p 32 (o)-p 23 (o))((7 + r)(2 7 -r + 6) 

+ 64J 2 ) 1 F 1 ( A_; 2 --;--) + 7(7' - e 2 )i^i («_; 1 - -;-- 

V 7 7/ V 7 7 

x ((p 23 (0) - p 32 (0))(r - e) + 8z J(p 22 (0) - 033(0)))), 
£ = 1.F1 («-; l - -; --) (-2 7 e(7 2 - e 2 )i^i («+; 1 + -\ -- 

V 7 7/ V 7 7 
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(A3) 

where e = \/r 2 — 64 J 2 , k± = -^ — 5 ,A± = -^ — -,»7± = 2^ and ii*i(a;&;z) is the Kummer confluent hypergeo- 
metric function [39] . These solutions are obtained by using Mathematica program. 
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